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S: 

I . We reported that the pair-pair correlation function of the two-dimensional 

O ' t — J model does not have long-range d-wave superconducting correlations in 

^ • the interesting parameter range of J/t < 0.5. The power-Lanczos method is 

> _ 
^^ I used under the assumption of monotonic behavior. This assumption has been 

T-H ! 

well checked in the two-dimensional t — J and attractive Hubbard model. Here 

^H ■ 
^- 

0^ . we re-examine this criterion of monotonic behavior of the pairing correlation 



function for the one-dimensional and two-leg t — J ladder where other accurate 



^ ■ numerical results are available. The method seems to be working well. 



PACS number: 74.20.-z, 71.27.+a, 74.25.Dw 



r> ' It has been shown that several essential physical properties of the high temperature 



superconducting (HTSC) cuprates can be described by the two-dimensional (2D) t — J 
model |ll|,0|, such as the single hole dispersion relation 0, non- Fermi liquid behavior [Q, 
and phase separation 0. One of the critical issues is whether the model gives large enough 
pairing correlation to quantitatively explain the high Tc. We reported part of the numerical 



results on this issue in Ref. which concludes that the 2D t — J is not superconducting 
in the physical interested region of electron density Ue and coupling constant J/t from the 
studies on the pair-pair correlation function and the two- hole binding energy [0. 

The criterion to evaluate the pairing correlation of the 2D t — J model used in Ref. 
I^J^ is to apply the power-Lanczos(PL) method [^ on the variational trial wave functions. 
We denote the initial trial wave function as | PLO) and the one improved by one Lanczos 
iteration as | PLl). We calculate the change of the pairing correlations of | PLl) and 
I PLO) for various variational parameters. If we assume the behavior of the long-range 
pairing correlations which is monotonic through the Lanczos or power iteration, then the 
one whose pairing correlation is least changed gives the correct results of the ground state. 
It is extremely difficult to prove this criterion by the present computing power in 2D, we 
could only provide more examples to support it. In a previous paper we have checked the 
validity of the criterion of the 2D attractive Hubbard model [^. Here we will show the 
results of the one-dimensional (ID) t — J model and 2-leg t — J ladders. 

The t — J model is: 

H = -t Y. i^tcja + h.c.) + J ^ (Si ■ Sj - -niUj), (1) 

<i,j>(T <i,j> 

where c]^ = q^(l — Ui^^), and (z, j) denotes the nearest neighbors i and j. 

For the ID case, we study the pairing correlation function in the spin-gap region of the 



phase diagram, i.e., large J(> 2) and low electron density [|l^,|TT|]. Here we use the correlated 



spin-singlet pair(CSP) wave function as | PLO). The CSP wave function is defined as 



\CSP) = PdU 



i>j 
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[E h-%f^/' |0) , (2) 



n=l 



where h = 2t/J. N^. is the total number of electrons, Pd is the projection operator that 



forbids two particles occupying the same site and the operator bl^ = Z^i Ci|c|^„| — c||cj_,_„|. 
The variational parameter u control the long-range correlations. In Fig. 1 we show the 
pairing correlation function of the | PLO) and | PLl) with J/t = 2.8. In order to check the 
ansatz of monotonic behavior, we choose a small lattice of 24 sites with 4 electrons which 
can be exactly diagonalized by Lanczos iteration. As it clearly show that the one with 
z/ = 0.1 has the stablest long-range correlation which is also very close to the value of the 
exact ground state. For smaller u the long-range correlation decrease while the larger ones 
increase. The results for this ID case are consistent with our ansatz. 

For the ladders discussed here, the leg-leg and rung-rung t and J are assumed to be 
identical. In the variational level of the 2 x L^ t — J ladder, the trial wave function used 



here is |12,13]: 



L 



HA, u)) = \l{-sin{j{x. - x,)Y ■ P, \{{u, + Sfccl^^cl, J | 0) (3) 



with Vk/uk = Afc/(efc + J el + Af,), A^ = A{coskx — 2cosky) and e/c = {2coskx + cosky) — /i. 
A is the d-wave superconducting order parameter and /x is the chemical potential. The 
operator Pd enforces the constraint of no double occupancy. We take t = 1 in this paper. 
In the 2-leg ladder case, the rung-rung pairing correlation P{R) = (A|Aj+/j) is measured, 
where A^ = (9,2^9, i;| - Cj-,2;TCi,i;i). 

For the ID and ladder cases, because the correlation functions decay exponentially or in 
power law with respect to the distance R, rather than being flat plateaus in two dimensions, 
it is difficult to define the "stable" correlation of the trial wave functions with different 
parameters A and z/. Thus we define the deviation D{A,v) to quantify the difference 
between the | PLO) and | PLl) wave functions. The wave function with minimal D is the 



one having stable correlation. It is defined as 

D = ^ E{(^i(^) - Po{R))/Po{R)r (4) 

where Po{R) and Pi{R) are the rung-rung pairing correlations, respectively. Note that this 
is not the unique way to define the deviation. We tried some other reasonable definitions 
of deviation and the results are consistent. Fig. 2 shows the wave function with minimal 
deviation of 12 electrons (61 and 6|) in the 2x8 ladder and J/t = 1. The wave function 
with minimal deviation is $(A = 0.12, z/ = 0.6). Fig.3 shows Po{R) and -Pi(-R) of the trial 
wave functions and the exact results of the system. The "sandwich" behavior of the long 
range part (R=3, 4) is clear. 

Since 2 x 8 is too small to see the long rang behavior of the pairing correlation, we did 

the similar analysis on 2 x 30 system, which the DMRG results are available for comparison 

g. We tried 24 t and 24 | electrons in 2 x 30 ladder with J/t = 1.0. The best one 



(minimal deviation) is A = 0.16 and z/ = 0.3, while the trial wave function with lowest 
energy is A = 0.24 and u = 0.1, whose pairing correlation is evidently overestimated. The 
pairing correlations obtained with three different sets of parameters are compared with the 
DMRG data from [jl4| in Fig. 4. It can be seen that the slopes in the log-log plot, that is. 



the exponents of the power law decay, of the trial wave function (-0.835) with minimal D 
and the DMRG data (-0.840) are very close. The differences of amplitudes and slopes of the 
long tails between these two may come from the different boundary conditions of DMRG 
(open boundary condition, OBC) and power-Lanczos method (periodic boundary condition, 
PBC). In fact, we have seen similar difference between OBC and PBC by solving the 2x8 
ladder with 12 electrons and J/t = 1 exactly, which is shown in the inset of Fig. 4. 



In summary, from the study of ID and 2- leg ladder t — J models, there are more evidences 
to support the "least deviation" criterion used to determine the superconductivity of the 
2D t — J model. It seems that the pure 2D t — J model doesn't give enough amplitude of 
pairing correlation in the physically interested parameter region. Since there are a number 
of evidences that the 2D t-J model is a fairly good model for HTSC, it seems the model 
should not be excluded just due to the negative result with respect to the d-wave pairing 
correlation. It is possible that we can keep the good properties of the 2D model and enhance 
the superconductivity by some reasonable modification, for example, adding the interlayer 



hopping term [16 
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Figure Caption: 

Fig.l: Pairing correlation function P{R) for various trial wave function(PLO) and their PLl 

results, together with the exact data(GS) for J/t = 2.8 of 4 electrons in a 24-site chain. 

Fig.2: Deviation D{A, z/) for a 2 x 8 t-J ladder, J/t = 1. 

Fig. 3: Pairing correlation function P{R) for trial wave functions with A = 0.12 and various 

u and their PLl results, together with the exact data for J/t = 1 of 12 electrons in a 2 x 8 

ladder. 

Fig. 4: Similar to Fig.2, but for 48 electrons in a 2 x 30 ladder. The dot lines shows the power 

law fitting. The inset shows the 2x8 exact results for the different boundary conditions. 
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